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^ < Abstract 

3 ■ 

Let u[x,t) be a (possibly weak) solution of the Navier - Stokes equations on all 

of R , or on the torus R /Z . The energy spectrum of «(-,?) is the spherical 
integral 

E{K.t)= I \u(k,t)\^dS{k), 0<k:<oo, 

■'\k\=K 

,J^ , or alternatively, a suitable approximate sum. An argument involking scale invari- 

CIh' ance and dimensional analysis given by Kolmogorov [7, ^ and Obukhov 1141 

predicts that large Reynolds number solutions of the Navier - Stokes equations 
in three dimensions should obey 



E(Kj)r^Ci)E^I'^K-^l'^ 



over an inertial range K\<k< K2, at least in an average sense. We give a global 
estimate on weak solutions in the norm ||=^5j:"('if)||~ which gives bounds on 

Cn ' a solution's ability to satisfy the Kolmogorov law. A subsequent result is for 

1^ , rigorous upper and lower bounds on the inertial range, and an upper bound on the 

"^ j ■ time of validity of the Kolmogorov spectral regime. © 2000 Wiley Periodicals, 

in ■ ^"'^• 

o 
oo 

O 



1 Introduction 



An important issue in the study of solutions of the Navier - Stokes equations in 

, ^ ■ the large is the principle governing the distribution of energy in Fourier space. The 

^ theory of Kolmogorov ITJUUll and Obukhov [14] plays a central role, predicting 

H power law decay behavior of the Fourier space energy density for solutions which 

- - exhibit fully developed turbulence. In outline, a basic prediction is that energy 

spectral functions E{K,t), or possibly its average over a statistical ensemble, is 

expected to satisfy 

(1.1) E{K,t)c^CQe^l'^K-^l^ 

over an inertial range of wavenumbers k" G [k"i, K2], where Co is a dimensionless 
constant, e is a parameter interpreted physically as the energy transfer rate per unit 
volume, and the exponents are determined by dimensional analysis [14J[6J. This 
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famous statement has been very influential in the field, and considerable experi- 
mental and numerical evidence has been gathered to support it. Despite its success, 
there have been relatively few rigorous mathematical results on the analysis of solu- 
tions of the Navier - Stoke equations, with or without bulk inhomogeneous forces, 
which have addressed the question as to whether solutions exhibit spectral behav- 
ior as described by (ll.ll ). Among those papers which do address certain aspects 
of these questions, we cite in particular two sources. Firstly, the book by C. Do- 
ering & J. Gibbon fSl reviews the Kolmogorov - Obukhov theory, and discusses 
the compatibility of spectral aspects of solutions with the L^ regularity theory for 
the Navier — Stokes equations. Secondly, S. Kuksin ifTTI proves that solutions 
to the nonlinear Schrodinger equation with added dissipation and with stochastic 
forces exhibit spectral behavior over an inertial range, with some positive exponent 
(which is not known explicitly). This latter work serves as an important mathemat- 
ical model of generation of spectral behavior of solutions under stochastic forcing, 
despite the basic difference in the equations that are addressed. 

In this paper we give a new global estimate in the norm ||=^5xM('>0lk°° o"^ 
weak solutions of the Navier - Stokes equations which have reasonably smooth 
initial data and which are possibly subject to reasonably smooth inhomogeneous 
forces. This estimate has implications on the energy spectral function for such 
solutions, and in particular in the case that there is no inhomogeneous force, we 
show that weak solutions of the initial value problem have spectral energy function 
which for all K G M+ satisfies 

(1.2) E{K,t)<AnR\, 
and time averages which satisfy 

(1.3) - / E{K,t)dt<^^, 

again for all k, where v is the coefficient of viscosity. In the case that a bounded 
inhomogeneous force is applied to the solution of the initial value problem, we find 
similarly that 

(1.4) E{K,t)<ATiR\{t), 
and furthermore 

(1.5) - / E{K,t)dt< 2V ; 



T Jo ' ' ' ~ VK^T 
In the situation of forcing being given by a stationary process, it is to be expected 
that the quantity R\{T)/T has a limit /?2 for large T , giving a constant upper bound 
for the time average of E{k,T). Since these estimates give rigorous bounds on 
E{K,t) with a faster rate of decay in wavenumber K than (11.11 ). this result presents 
a conundrum. Either it is the case that solutions which exhibit large scale spectral 
behavior as in dl.ll ) are not smooth, and in particular do not arise from the initial 
value problem with reasonably smooth initial data. Or else the bounds (1 1.21 ) (1 1.31 ) 
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(and (II .41 1 (1 1.51 ) respectively, in the case with inhomogeneous forces) give restric- 
tions on the spectral behavior of solutions, and in particular an upper bound on the 
value of the parameter e, a restriction on the extent of the inertial range [k"i, K2], 
and in the case of (1 1.21) (II .31 ). an upper bound on the time interval [0, Tq] over which 
spectral behavior may occur for a solution. 

There is a well developed literature on the energy transfer rate e, and other as- 
pects of the Kolmogorov - Obukhov theory, based on physical assumptions on the 
character of the fluid motion. These assumptions are for flows exhibiting fully de- 
veloped turbulence, and are described in Obukhov [14], for example. They include 
the hypothesis that the flow is in a stochastically steady state, the energy transfer 
rate is of a certain form and exhibits a particular scale invariance, and that the en- 
ergy spectral function is negligible for wave numbers higher than a cutoff Ky. Un- 
der these assumptions, the cutoff Ky is determined (it is known as the Kolmogorov 
scale) and the energy transfer rate e is identified with the energy dissipation rate 

(1.6) e,:=^^[\-E{K,t)dK. 

Using the latter statement, that £ = £1 , it is possible to find a better estimate of 
the behavior of £ with respect to Reynolds' number than ours in this paper, as 
for example in C. Doering & C. Foias [4J. The difference between this body of 
work and our analysis is that we make no physical assumptions on solutions of the 
Navier - Stokes equations, deducing our conclusions purely from known results 
about such flows. It is worthwhile to point out as well that the upper and lower 
bounds to the inertial range in our work are conclusions, as compared to prior 
work in which the upper bound Ky on the inertial range is an assumption of the 
theory, and no explicit lower bound is given. 

In Section [2] we give a statement and the proofs of our estimates on the Fourier 
transform of weak solutions of the Navier - Stokes equations, posed either on all 
of X G M^ or else for x G T^. Since there is no known uniqueness result, one cannot 
speak of the solution map for Navier - Stokes flow, and we emphasize that this es- 
timate is valid for any weak solution that satisfies the energy inequality. In section 
[3] we interpret these estimates in the context of the spectral energy function, and we 
analyse the constraints on spectral behavior of solutions mentioned above, giving 
specific and dimensionally appropriate estimates for the endpoints of the inertial 
range K'l, K2. In the case of no inhomogeneous forces, we give an upper bound 
To on the time of validity of the spectral regime. The bounds on K^i and K2 are 
also valid in the probabilistic setting, for statistical ensembles of solutions. That is, 
suppose that one is given an ergodic probability measure (P,^) on the space of 
divergence free vector fields which is invariant under some choice of definition of 
Navier - Stokes flow. As long as the inhomogeneous force and the support of the 
invariant probability measure are contained in the closure of the set of reasonably 
smooth divergence free vector fields, then our constraints on the spectral behav- 
ior of solutions apply. The final section gives a comparison of our constraints on 
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{ki,K2,Tq) to the Kolmogorov length and time-scales of the classical theory, and a 
discussion of the dimensionless parameter ry := k:2/k'i as an indicator of spectral 
behavior of solutions. 



2 Estimates on the Fourier transform in L°° 

The incompressible Navier - Stokes equations in their usual form are written 
for the velocity field u{x,t) of a fluid, its pressure p{x,t), and a divergence-free 
force f{x,t), 

dtU + {u-V)u = -V/7 + vAm + / 

(2.1) V-M = 0, 

m(-,0) =mo(-) , 

where we consider spatial domains either all of Euclidian space ;c G M^, or else 
the compact and boundary less torus a; G T^ := M^/F, where F C R^ is a lattice of 
full rank. Denote by D either of the above spatial domains. The time domain is 
< f < +00, and the inhomogeneous force function / is assumed to be divergence- 
free and to satisfy / G LJ^,^([0,+oo);//-i(D) nL^(D)). A 'Leray' weak solution to 
(12.11 ) on D X [0, +oo) satisfies the three conditions. 

(1) Integrability conditions: For any T > the vector function {u,p) lies in 
the following function spaces 

(2.2) MGL~([0,r);L2(D))nL2([0,r);//^(D)) , 

(2.3) peL'l\Dx[0,T)), 

(2) Weak solution of the equation: the pair {u,p) is a distributional solution of 
( 12.11 ). and furthermore lim,^o+ u{-,t) = uq{-) exists in the strong L^ sense, 

(3) Energy inequality: the energy inequality is satisfied 

(2.4) \ / \u{x,t)\^dx + v I I \Wu{x,s)\^ dxds 

Jd Jo Jd 

— u{x,s)- f{x,s)dxds < J / |mo(x)| dx 

Jo Jd Jd 



for all < f < +00. The inequality (12.41 ) is an identity for solutions which are 
regular. It is well known that weak solutions exist globally in time, either when 
/ = 0, a result due to Leray fTT.TTl, or when / is nonzero. The question of their 
uniqueness and regularity remains open. 

Many facts are known about weak solutions, including that for any T > the 
interpolation inequalities hold; u G U{[0,T);LP{D)) for all 3/p + 2/s = 3/2, for 
2 < /? < 6. That the L^/^ estimate for the pressure in ( |2.3b is sufficient is due to 
ifTSll . Considering a weak solution as a curve in L^{D) defined over t G M+, the 
following proposition holds. 
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Proposition 2.1. A weak solution is a mapping [0,r) h^ L^{D) satisfying the con- 
tinuity properties 
(2.5) 
ui.)GL-i[0,Ty,LHD))nC{[0,Ty,L\D)^,^,,i^p^l^gy)nCi[0,Ty,H-'{D)) 

for any 5 > 0. Furthermore, as a curve in Sobolev space, 

^GL\[0,Ty,H-nD)). 

A clear exposition which includes these basic facts is the lecture notes of J.-Y. 
Chemm Q. 

Being a curve in L}{D), the Fourier transform of a weak solution makes sense, 
andM(-,0 = ^M(-,0isagainacurveinL~([0,r);L2)nC([0,r);L2^^^j^p^j^gy). 

We will make use of a dimensionally adapted Fourier transform ^, namely 

(2.6) u{k) = {^u){k) := -^ I e-'^''u{x)dx , 

V^l^ Jd 

where /: G M^ when the spatial domain is D = M?, and we set V = (27r)^ in stan- 
dard units of volume in M^. In this setting, the norm of the Fourier transform 
is given by \\u\\^ = J^3 \u{^)\'^d^. When the domain is D = T^ = R-^/F, we 
take k £ V the lattice dual to F, we set V = |F| := vol(M^/F), and we define 
II "IP •= HkeV l"(^)P|r'|- With this choice, the Plancherel identity reads 

„ „9 V „ „9 

(2.7) ll.ll =(^ll«ll • 

With respect to the normalization, the function u{x,t) has units of velocity L/T, 
and its Fourier transform is such that \ii{k,t)\-^ has units of Fourier space energy 
density {L/TfL^. 

2.1 An estimate on ^dxu{-,t) on the torus 

Focus the discussion on the case of the spatial domain D = T^. Then any choice 
of initial data uq{x) G L^(T^) has uniformly bounded Fourier coefficients, indeed 
|"o(^)| < II Wo II- Furthermore, since the complex exponential e^'^'^ is a perfectly 
good element of {H^^'^)* which, being tested against u{-,t) gives the Fourier co- 
efficients, we also have the result 

Proposition 2.2. For each /: G F' the Fourier coefficient u{k,t) G C^ is a Lipschitz 
function oft£ M+. 

It is again made clear by this that the problem of singularity formation is not 
that u{k,t) becomes unbounded, but rather that H^ mass, including possibly L^ 
mass, is propagated to infinity in ^-space in finite time. 

A (future) invariant set A is one such that mq G ^ implies for all f > 0, u{t) G A 
as well. When / = the energy inequality (12.41) can be viewed as implying that the 
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ball Br(0) C L}{D) is an invariant set for weak solutions satisfying ||mo||l2|-£,) < R. 
In similar terms, we give another invariant set for weak solutions. Define the set 

(2.8) A«, ■={u : VA:Gr',|/t||M(/t)| </?i} , 
and as above let Bfi(O) denote the ball of radius R in L^(T-'). 
Theorem 2.3. In the case that / = 0, whenever 

(2.9) -= < vRi 



then the set Ar^ nBR{0) is invariant for weak solutions of the Navier — Stokes 
equations (12. lb . Thus, if the initial data satisfies uq € Ar, 08^(0), for allO < t < 
+00 the Fourier coefficients of any Leray weak solution emanating from this data 
satisfy 

(2.10) sup \u{k,t)\ < -i- , Mk£T' . 
0<f<+<» I'^l 

This result appears in the paper m in a slightly different form. For nonzero 
/ the ball Br{Q) C L^{D) is not necessarily invariant. However given uq G Br(0) 
and our hypothesis that / G L^^([0,+oo);//^i(D) r\L?{D)), there is always a non- 
decreasing function R{T) > R such that for all T > 0, u{-,T) G B^j~^{Q). Indeed, 
suppose that a Galilean frame is chosen and the pressure p is suitably normalized 
so that Jj)u{x,T)dx = = fj)f{xj)dx. Let F'^{T) := J^ ||/||^_i dt. Then by 
standard interpolation one has that 

(2.11) \\u{.,T)\\l, + vJ^ \\Vu{;s)\\l,ds<R\T). 

The function R^{T) is an upper bound for the LHS of the energy inequality, for 
which there is an estimate R^{T) < R^ + ^F^{T). In case of a bounded inhomoge- 
neous force /(•,?) G L°°(R+;//^^ nL^), there is an upper estimate F^{T) < CT, so 
that/?^(r) exhibits (not more than) linear growth in T. 

Theorem 2.4. In the case of nonzero f{x,t), let R{t) be an a priori upper bound 
for ||M(-,f)||^2, for example as in (12.111 ). Suppose that R\_{t) is a nondecreasing 
function such that for all {k,t) we have 

(2.12) __ + ^_<v/?i(0, 

then the set {{u,t) : 0<t, u(-,t) GAjf^/AnBi^tf^iO)} is invariant forweaksolutions 
of the equations (12.11) . That is, if the initial data satisfies uq G A^j/oj nBs(o)(0) t^^^^ 
for alio <t < +00 the Fourier coefficients of any weak solution emanating from uq 
and being subject to the force f will obey the estimate 

(2.13) |^(^,,)|<^. 
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It is a common situation for the inhomogeneous force to have properties of 
recurrence, such as if it were time-periodic, or if given by a statistical process which 
is stationary with respect to time. For bounded / as above, the estimate F^{T) < 
CT holds. Furthermore, one is interested in those solutions which are themselves 
statistically stationary. For these solutions it will be the case that the force adds to 
the total energy at the same rate as the dissipation depletes it. Indeed, one assumes 
that there is a constant R such that ||M(-,r)||^2 < R and v/q ||Vm(-,5)||^2<5?'5' < 

2 

R T. Therefore, in the statistically stationary case we expect that the upper bound 
R{t) in the hypotheses of Theorem 12.41 to be given by the constant /?, which in 
particular gives a uniform bound in time. In this situation, the constant 7?i is also 
time independent. 

Proof of Theorems \Z3\ and \Z4\ For each keV the vector u{k) € C^ ^ C^, where 
C| = {w G C^ : w -Lk = 0} is specified by the divergence-free condition. Because 
{u,p) is a distributional solution, the Fourier coefficients satisfy 

(2.14) dtii{k,t) = -v\kf-u{k,t)-iUkk-—=Y^u{k-ki)®u{ki)+f{k,t), 



^i, 



X{u)k 



at least in the weak sense, after testing with a smooth cutoff function <p(f) G C°°. 
We use the notation that X(m)^ is the k'^ component of the vector field represented 
by the RHS. The convolution has introduced the factor \/\/V . The operator 11^ : 
C^ -^ Cl is given by 

k 

\k\ 



U,{z)=z-{z-k)- 



The Leray projector onto the divergence-free distributional vector fields, consid- 
ered in Fourier space coordinates, is the direct sum of the Ilk- 

The radial component of the vector field X{u)k in C\ C C^ is expressed by 
xe.{u{k) •X{u)ji)/\u{k)\. Consider first Theorem 12.31 which is the case that / = 0. 
Suppose that |/:||m(^)| = R\ for some k, that is, the solution is on the boundary of 
the region A« I . Since \Ylkk-Y,ki^{k — k])®u{k\)\ < |A;|||m||^2. an estimate of the 
radial component ofX{u)k is that 

re{a(]c)-X{u)k) = -v\k\^\u{k)\^ + -^im{a(k)knkY,u{k-ki)0u{ki)) 

<-vRJ + ^\\u{;t)\\l,R,. 



When u{-,t) £6^(0) and R^ < vRi \/V the RHS is negative, implying that integral 
curves {i{k,t) cannot exit the region. Thus the ball of radius Ri/\k\ in C^ is a 
trapping set, or a future invariant set, for the vector field X{u)k. 
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In case of the presence of a force /, suppose again that |^||M(A:,f)| = Ri{t) for 
some {k,t). The radial component of ^(m)^ satisfies 

re{a(k)-X{u)k) = -v\kf \u{k)\^ + -^im [0(1) knkY,ii{k-ki)(g)u{ki)) 

vV I, 



(2.15) +xe{u{k)-f{k,t)) 

<-vR\ + ^\\u{.,t)\\l,R, + \f{k,t)\^^. 

Furthermore, the energy 1 1 w ( • , ? ) 1 1 ^2 is bounded by R^{t). As long as the radial com- 
ponent ofX{u)k (namely the quantity in (12.151 ) normalized by the length \u{k)\ = 
R\ {t)/\k\) is bounded above by the growth rate of the ball itself, namely by Ri/\k\, 
then solution curves (u{-,t),t) do not exit the set {(«(•),?) : < f , u{-,t) G 
A/jj(f) nS/j(j)(0)}. In particular this happens for nondecreasing R\{t) whenever 
{^)-'R\t) + \f{k,t)\/\k\<vR,{t). 

Therefore when the initial data uq G ^Ri{o) and the force / satisfies (12. 12I ). then 
the solution satisfies u{x,t) G A^|(f) iox R\{t) finite, for all positive times t. D 

It is natural to ask what constraints are imposed on the data uq by the condition 
(12.91 ). Given smooth initial data mq and a force / satisfying \f{k,t)\ < F2\k\, the 
constant vRi can always be chosen so as to satisfy ( 12.91 ) (the case / = 0), or if the 
force / 7^ 0, the function Ri{t) can be chosen to be nondecreasing and to satisfy 
(12.121) . Thus the hypotheses to this theorem encompass any reasonable smooth 
class of initial data and inhomogeneous forcing terms. We note that the constant 
Ri scales dimensionally in terms of the units l}!'^ jT . 

Under changes of scale, the quantity supySup^|/:||M(/:,f)| transforms like the 
BY-nonxi sup^ ||5v"(')0IIl1' and indeed the latter being finite implies the former. 
However as far as we know there is no known global bound on the BW-nonn of 
weak solutions to (12.11 ). A related inequality appears in [3|, which is a global upper 
bound on the L^-norm of the vorticity co := V x u, again uniformly in time. 

A corollary to this result gives a stronger estimate for time integrals of the 
Fourier coefficients of weak solutions. 

Theorem 2.5. Suppose that a weak solution u{x,t) satisfies (12.111 ). with its initial 
conditions satisfying uo G Ar,. Then time integrals of the Fourier coefficients obey 
the stronger estimate 



rT 

(2.16) / \u{k,t)\'^ dt < 

Jo 

The constant is given by 



Rl{T) 



R2{T) = \ (r,{T) + ^2R\{0)+Rl{T)^ 
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vVV vv V>/() 



Assuming that supi^^^3Fi{k,T) := Foo{T) < +00, ;/je constant RiiT) will be inde- 
pendent ofk. 

When the force / = 0, the constants R and R\ can be taken independent of T , 
implying that R2 is also constant in time, and the estimate (12.161) holds uniformly 
over < r < +00. For nonzero forces which are L°° with respect to time, there 
is an upper bound Fi{k,T) r^ \/Y, and R2{T) will grow at most linearly in time 
for large T. In the situation in which the forcing is given by a stationary process, 
it is expected (but not proven at this point in time) that for typical solutions, the 
quantity R\{T)/T will have a limit /?2 for large time T , representing a balance 
between energy input and dissipation. Notice that R2 scales dimensionally in terms 

ofLV2/r. 

Proof. Because the field u{-,t) is divergence-free, k •u{k, •) = 0, implying the vector 
identity k ■ u{k — k\)® u{k\) = u{k — k\)-k\® u{k\). The absolute value of u can 
be estimated from (12.141 ). 

\dt\u{k,t)\^ + v\k\^\u{k,t)\^ = im^={u{k,t)-nkY^u{k-ki,t)-ki®u{ki,t)) 



+u{k,t) ■ f{k,t) , 

which is valid for each k in the sense of weak solutions in time. When integrated 
over the time interval [0, T] it gives 

fT 

(2.17) v|A:|M \u{k,t)\''-dt = ^\uQ{k)\'^-\\u{k,T)\^ 

I rT 

+im— = / (ii{k,t) ■Ilk'Vi2{k — ki,t)-ki0ii{ki,t))dt 
VVJo t' 



+ / u{k,t)-f{k,t)dt 
Jo 



Multiplying this identity by \k\^, the terms of the RHS are then bounded as follows: 

'j\k\^\uo{k)\^ < ^/??(0) , 



\k\^\ f u{k,t)-f{k,t)dt\<-^(v\k\^ f \u{k,t)\^dty'^{l \f{k,t)\^dty^^ , 

Jo yv JO JO 

and finally 

I rT 

1^1 ^^lim / (u{k,t) ■Tliiy\u{k — kY,t) ■ki(^ii{ki,t))dt\ 
VV' Jo "i^ 

<J—{v\k\'r\u{k,t)\^dtf^ sup \\u{;t)\k2{v r\Vu{;s)\^dsf\ 

vVV Jo o<t<T Jo 
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Define f{k,T) = v\k\^ Jq \u{k,t)\'^dt, then the identity ( 12.171 1 implies the inequal- 
ity for /(A;,r) 

p,8) ,^(„)<^ + (^ + ^);(.,r). 

' ' 2 V vVv Vv ^ 

where we have used that R{t) is nondecreasing. The quantity I{k,T), being non- 
negative, cannot exceed the largest positive root of the quadratic equation where 
equality is attained, giving the estimate ( I2.16I ). This estimate is uniform ink ^1? 
as long as sup^^^s Fi{k, T) = F„{T) < +oo. We note that 

keZ3\{o} II 
which appears in the discussion of energy estimate bounds ( I2.11I) . D 

2.2 The analogous estimate on M^ 

Suppose that ||m(-,0IIl2 ^^(0 (if there is no force, then/?(f) = R{0) suffices). 
The main difference in the case of D = M^ is that the functions u{k,t) are ele- 
ments of a Hilbert space, whose values at a particular Fourier space-time point 
{k,t) are not well defined. We work instead with filtered values of the vector 
field u{x,t). LetO^ke E?, and for 5 < \k\/{2y/3) define Xki^) a smooth cut- 
off function of the cube Q^ about A; G M^ of side length 25 (5 < 1 is acceptable 
for large \k\) which takes value ;^^ = 1 on a cube of half the sidelength. The point 
is that for ^ G supp(^i) then \k\/2 <\^\< 3\k\/2. Now define {XkiD)u){x,t) = 
^^^Xk{^)u{B,,t) = {Xk*u){x,t). Since Xk G H'" for all m, it and its translations 
are admissible test functions, the statement (12.51 ) implies that {Xk* u){x,t) is a 
Lipschitz function off for each x. Define ep{k,t) := {J\Xki^)u{^,t)\P d£,y^'^ for 
2 < p < +00, the conclusion is the following. 

Proposition 2.6. The function ep{k,t) is a Lipschitz function oft £ M+. 

We quantify the Fourier behavior of the force / in similar terms. Consider the 

functionz.(D)/(x,0 = ^-'z^(^)/(^,Oandlet/^(^,0:=supo<,<,(/(|z.(^)/(^,OIVI^I^)^§)'/^- 
Recalling that/ E Lr,,([0, +oo);//-i nL2(D)) we have /2(^,0 < supo<,<, ||/(^,s)||ri-, . 
However the fact that fp is finite is in general additional information about the reg- 
ularity of the forcing. 

Theorem 2.7. Let initial conditions uo(x) give rise to a weak solution u{x,t) which 
satisfies \\u{-,t)\\ii < R{t). Suppose that there exists a nondecreasing function 
R\ {t) such that for all 2< p < +oo and t E M+ 

(2.19) ^25fp^+f^^(^k,t)<Uyit), 
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where 5 < |/c|/2v3- Consider a solution to ilAi that initially satisfies sup2^ ., ^ , „ g^, (^ . 0) < 

Ri{0)/\k\. Then for all teW+ 

(2.20) |^,(^)^(^,,)|^^<^. 

Theorem 2.8. Suppose a weak solution u{x,t) satisfies (12.111 ). and furthermore ask 
that \Xk{^)uo{^)\L- < ^i(0)/|^|. Then for all T G M+, 

(2.21) f^\M{;t)\l^dt<^^, 
where the constant R2{T) is given by 

(2.22) R2{T) = \ (RsiT) + ^4/?2(0) +/?2(r)) , 

where 
(2.23) 

R5{T)=^-^+^-^, F^{T)= sup (rm^)M,t)\i.dty^\ 



V\/V VV k(=Mi\!n\^Jo 



The strategy of the proof of these two results is to give an analysis similar to 
that of Section l2?T] for a uniform bound on ep{k,t) with the correct behavior in the 
parameter k. The first lemma controls the behavior oie2{k,t), pointwise in t. 

Lemma 2.9. Suppose that R\ (t) is nondecreasing, and is such that for all t £ R+ 

(2.24) (^25)y^^+f2ik,t)<^Riit), 

where 5 < \k\/2V3. Ife2{k,0) <Ri{0)/\k\, then for all 0< t < +oo 

(2.25) e2{k,t) < ^ . 



Proof. The quantity e2{k,t) satisfies the identity 

{2.26^j^el{k,t) = ^drJ\m^d^ 

= re|(MM7))(-v|^P(Z.(^)«(^,0) 

-iU^)^n^Ja{^-^,,t).^,<g)u{^,,t)d^,)d^ 



12 A. BIRYUK & W. CRAIG 

The first term of the RHS is negative, bounded above by 

-rev J {Xki^)u{^,tm\\U^)ui^,t))di = -v J \^\'\M^)u{^,t)\'d^ 

where we recall that \^\ > {\k\ — \/35) > \k\/2 holds for ^ £ supp (Xk)- The sec- 
ond term of the RHS of (12.261 ) is bounded with two applications of the Cauchy - 
Schwartz inequality; 

(2 



<j^Uku\yUkn^^-Jui^-^ut)^ii{^i,t)d^iy 



yV f 

< j^^UkuWi^ KMl^ \J u{^-^i,t)0u{^i,t)d^^^ , 

where we have used the property of incompressibility that u{^ — ^i) ■ ^i = t,-u{t,— 
t,\). Furthermore on the support of ^^, |(^ | < 3|^|/2 therefore 

\MkU < ^(25)^/2 , \ju{^-^ut)m(^,,t)d^y I,. < pii^, < ^2£}^ 

The third term of the RHS of (12.261) is not present without a force. When there is a 
force, it admits an upper bound 

\ve jjMM^)u^)f{^.t)di\ < mmy \m-'xJi;t)\\L^ 

<^IIZHlL^II/(-,Ob-. ■ 
An estimate of the RHS is therefore 

RHS < -^\k\hj{k,t) + IM^-L R2^t) \k\e2{k,t) + ^f2{k,t) \k\e2{k,t) . 

This is the situation from which the proof of Theorem l2.4l proceeds. Consider the 
set B«, = {e : e < {Ri/\k\)}, and suppose that the inequality holds 

(2.27) ^^ R\t) +f2{k,t) < -Ri{t) . 

yJV 6 



When e = ^2 is on the boundary of B«, , that is when 62 = Ri {t)/\k\, then 

V 3f25)^/^ 1 3 

RHS < --\k\^el{k,t) + ^ ^' -= R^t) \k\e2{k,t) + -f2{k,t) \k\e2{k,t) 



< (-^/?i + (25)^/'-^ R\t) +f2ik,t)) ^-R, < 
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That is e2{k,t) < 0, and thus Br, is an attracting set for e2{k,t). If initially ^2(^,0) < 
Ri{0)/\k\, then for all t € M+, e2ik,t) < Ri{t)/\k\. This proves the lemma. D 

Lemma 2.10. Given k G M?, suppose that for some 2 < p < +00 there is a nonin- 
creasing function R[{t) which satisfies 

(2.28) (^25)yP^+fp{k,t)<^R,it), 

for some 5 < |^|/2v3- If a solution to (12.11 ) initially satisfies ep{k,0) < R\ {0)/\k\, 
then for all t G M+ 

(2.29) ep{k,t)<^. 



Proof. The principle is to show that the local U' norms of u{^,t) are bounded, 
using the same strategy as the proof of Lemma [2!9l Since ep{k,t) is Lipschitz 
continuous for each ^ G M^, one calculates 



{230)j^ePik,t) = d,J\Xka\"d^ 



= re J p\MrHU^)K^,t))d,m^)ui^,t))d^ 

= Te J p\m'-'mm^{-y\^\\u^)H^,t)) 

+rc J p\M'-'lMm^)M^)fi^,t)d^ . 
The first term of the RHS of ( 12.301 ) is negative, 

-pvJ\^\^\Xku\'d^<-pv^-^ePik,t). 

Using the assumptions of the lemma and the Holder inequality, the second term 
has an estimate 



im I p\XkurHZk{^)u{^,t)) [M^)^n^^.ju{^-^ut)0u{^,,t)d^,)di 



p[j\Xku\'d^)''' '{JimXkm'd^) ''^\Ju{^-^ut)®u{^ut)d^ 



{2ny 



<\k\eP-'^-l{26)yP^RHt) 
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The third term of the RHS of (l230l) is bounded by 

An estimate of the RHS of (12.301) is thus 

Consider again the set B^i = {^ : < e </?i/|^|}. When e = e^ is on the boundary, 
that is when Cp =R\{t)/\k\, then 

V < 3(25)3/^ 1 2, -^r^ 3 /?r^ 






4|A:|P-2 2 VV \k\P-^ 2-"'\k\P- 



Supposing that ( I2.27I ) holds, the RHS is negative for e = Cp on the boundary, and 
the set Bsj is attracting for the quantity ep{k,t) for t G M+. D 



These are estimates which are uniform in the parameter p. We are now prepared 
to complete the proof of Theorem |2.7[ indeed we note that limp^+ooep{k,t) = 
IXkuli"^- The quantity ep{k,t) is given a uniform upper bound in Lemma l2. lOl under 
the stated hypotheses, and hence the theorem follows. 



Proof of Theorem 1251 Start with the identity in (12.301 ) for ep{k,t), which we read 
as 

(2.31) d4 = ^e f /'-''> d,ePp . 

Because of the support properties of the cutoff functions Xk{^), there is the com- 
parison \k\/2 < \^\ < 3|^|/2 on the support of Xk> thus one has upper and lower 
bounds 

(2.32) 
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We therefore can rewrite the RHS of (12.311 ) as 
(2.33) 2ef/'-'^(^-vl\^\^\U^)u{^,t)\Pd^ 



15 



+ 



{IkY 



im \Ui)u{^,t)r'M^)u{^,t 



xU^)(n^Ju{^-^,,t).^y^u{^,,t)d^y)d^ 



:=-/i+/2+/3 
The first of the three terms of the RHS is 



/i =2v 






^\'\U^)u{^,t)\Pd^ 



2/p 



which is well-defined because of ( I2.32K is positive, and through a lower bound will 
give us the result of the theorem. The second term of the RHS is 



n\M^)ui^,twd^ 



-1+2/ p 



for which one uses the Holder inequality (with {p — 2)/p + \/p + \/p = 1) to 
obtain an upper bound; 



2VV 



i/2i<^(/iz.«i^^^)^''"''/^(/iz.«r^^)^/^(/i^.r^^ 



-l+2/p 



X n^ju{^-^i,t)-^^®u{^^,t)d^i\^^(^j\Xku\Pd^ 
The third term of the RHS is 

h = 2v^ j \U^)u{£,,t)rHk{^)u{£,Mk{^)M.t)d£, 

/ f \-l+2/p 

x[J \Xku\''d^j 
which satisfies an estimate of similar form, namely 

=2{Jm\Pd^f''{J\xj\''dif''. 



-i+2/p 
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Integrating (12.311 1 over the interval [0,T], 



dt 



< el{k,0)-el{k,T)+l \h{l)\ + \h{t)\ill . 

Because of the properties of Xk, we have \k\/2 < |i^ | < 3|^|/2 in the support of the 
integrand, and therefore 



|)t|2-4/P<(2 



jmu^)u{^,t)\'d^ 



1-2/ p 



which means that the LHS of (12.341) gives an upper bound for the quantity 

Cancelling terms, one uses Cauchy - Schwartz to estimate the two time integrals 
on the RHS of (l234l) . Using that 



we estimate the first time integral as follows: 



X[£\\u{;t)\\l4^u{;t)\\l,dty 

<^i2^8Yil'ilmmf^d. 



1/2 



1/2 



x(^(sup \\u{;t)\y){f\\\Vu{;t)\\l,dtf') . 
V V ^^ 0<t<T -fO J 

We have used the Plancherel identity and its constant, as well as the fact that 

/|Zi|P^^< (25)3. Thus 



2/?2(r) 



Vp,\'I^ 



\i,\dt<::^^{2'5y/p(f v{f\m''d^f''dt 

vVV ^Jq J 



Under similar considerations. 



^''i/3iJ?<2^'(|i;e^«i''j^)'/^jf)'^'(^^(|i^jrc/§)'/"^f 



2/P,y/' 
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Now multiply the inequality (12.341 ) by \k\^ and use the above estimates with the 
fact that \k\/2<\t,\< 3|A:|/2 on the support of Xu. 

£ vlkl'ijim'd^f'dt <\k\mo) 

(2.35) +2{{25)y/P^[f\\k\\f \x,u\Pd^f''dty^' 

From our hypotheses on the initial data we know that |^pep(0) < R]{0). Defining 
fpik,T) := J^v\k\^{J\Zku\Pd^f^''dt, the inequality (12351) states that 

(2.36) Ilik, T) < RJiO) + (2(25)3/^^ + ^:l>^\(k,T) , 

^ V\/V \/v / 

where we define Fp{T) := (/q {J \Xkf\'' d^)^'^ dty'^ . As we have argued before, 
this implies that Ip{k, t) cannot exceed the largest positive root /?2,p of the associate 
quadratic equation, resulting in the statement that 

lp{k,t)<R2,p{T) 

where the constant R2,p{T) is given by 

(2.37) R2,piT) = i (/?5,p(r) + ^4/?2(0)+4^(r)) , 
where in turn 

%.(r)=2(25)^//'^ + ^, F'-{T) = [j\j\Mm^,trd^)y''dt). 

The result of the theorem will follow by taking the limit of large p — > +o°, recov- 
ering the estimate on IXk^^lL-"- D 

3 Estimates of energy spectra 

The energy spectral function is the main concern of the present paper. For the 
problem ( 12.11) posed on D = M^ this is defined by the spherical integrals 



(3.1) E{K,t)= f \u{k,t)\^ dS{k) , 

J\k\ = K 

where < k: < +oo is the radial coordinate in Fourier transform variables. When 
considering the case of a periodic domain D = T^ the Fourier transform is defined 
over the dual lattice, and therefore to avoid questions of analytic number theory 
one defines the energy spectral function to be a sum over Fourier space annuli of 
given thickness a; 

(3.2) ^('^'^) = ^ I^ |M(^,OI'|r'| . 

K<\k\<K+a 
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The classical Sobolev space nomis of the function u can be defined in terms of the 
energy spectral function, via the Plancherel identity. Indeed in the case D = M.^ the 
L^ norm is given as 

(3.3) Ml. = ^fEiK)dK, 
and the //'' Sobolev norms are 

(3.4) ||„||^. = ^^y^ {K^ + \yEiK)dK. 

Analogous definitions hold for the case x G T^. 

3.1 Kolmogorov spectrum 

There is considerable lore and a large literature on the behavior of the spec- 
tral function, particularly for large Reynolds number flows, the most well known 
statement being due to Kolmogorov. The prediction depends upon a parameter e, 
which is interpreted as the average rate of energy transfer per unit volume. Assum- 
ing that a flow exhibiting fully developed and isotropic turbulence has a regime of 
wavenumbers over which E{k,-) depends only upon e and k:, Kolmogorov's fa- 
mous argument states that over an inertial range k: G [k'i , K2], 

(3.5) £(^c,•)~Co£^/^fc-^/^ 

for a universal constant Cq. His reasoning is through a dimensional analysis. The 
actual history of this prediction, which is well documented in [6] among other 
references, includes a number of statements of Kolmogorov as to the small scale 
structure of the fluctuations in a turbulent flow |l7l[8l|9l, and an interpretation of his 
results by Obukhov lfT4l in terms of the Fourier transform, as is stated in (13.51) I 
Some of the issues surrounding this statement are whether the Kolmogorov scaling 
law (13.51 ) should hold for an individual flow at every instant in time, whether it 
should hold on time average, or whether it is a statement for the average behavior 
for a statistical ensemble of flows with the probability measure for this ensemble 
being given by some natural invariant measure for solutions of the Navier - Stokes 
equations. The bounds given below have implications on the energy spectral func- 
tion in all of these cases. 

3.2 Bounds on energy spectra 

The estimates given in section [2] on the Fourier transform of solutions translate 
into estimates on the energy spectral function for such solutions. Bounds which 
are pointwise in time are given in the following theorem. 



In fact Obukhov formulated an integral version of this result, which he called the 'two-thirds 
law of energy distribution'. 
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Theorem 3.1. Suppose that / = and that the initial data satisfies uq £ A^i n 
Bj{{0), where R and R\ satisfy ( I2.9I ). Then for all K and all times t, 

(3.6) E{K,t)<4KRJ. 

In the case of non-zero forcing f, then there is a finite but possibly growing upper 
bound given by 

(3.7) E{K,t)<47i:Ri{t) 

Bounds which concern the time average of the energy spectral function are 
derived from Theorem [ 



Theorem 3.2. Again suppose that the initial data satisfies uq G A^, nB«(0), where 
RandRi satisfy (|2!9l). and the force f £ Lfg^{[0,+oo);H^\D)nL^{D)) is bounded 
as in (12. 191 ) (12. 23 1 ). Then for every T the energy spectral function satisfies 



T AnRliT) 1 



(3.8) i/ E{K,t)dt< 

I Jo 



VT K 



■2 



In particular, under the hypotheses of Theorems 12. 3[ I2.4[ and 12.51 the energy 
spectrum must decay with an upper bound of order ^(k:^^) for every T. In case 
that the solution is such that limsupy^_|_„/?2(^)/^ is finite, then the time average 
behavior in (13.81 ) has an upper bound which is uniform in T. In any case, this is 
evidently faster than the Kolmogorov power law (13.51 ) and thus merits a further 
discussion. 

Proof of Theorems 13. l\ and U?2\ In the case of spatially periodic solutions, the def- 
inition of the energy spectral function gives that 

EiK,t) = 1 £ \uik,t)\'\r'\ 

K<\k\<K+a 
K<\k\<K+a I I 

We have used that the lattice point density of F' is |r'|^^ The inequalities of 
Theorem 13.11 follow. In the case in which the spatial domain D = M^, the proof is 
similar. 

To prove Theorem l3.2[ consider first the case of D = M^, where 

I- f E{K,t)dt = [ (^ r\u{k,t)\^dt)ds{k) 

T Jo J\k\=K^T Jo J 

VvrK-4. 
This is the stated estimate. The periodic case is similar. D 
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3.3 Estimates on the inertial range 

The two theorems 13.11 and 13.21 have impUcations on the inertial range of a so- 
lution of ( 12.11 ). In particular the inequalities ( I3.6N3.7I ) give uniform upper bounds 
ioT E{K,t), while (13.81 ) estimates its time averages from above with a decay rate 
Ck~^. For direct comparison we define the ideaUzed Kolmogorov energy spectral 
function with parameter £ to be 

(3.9) EK{K)=Coe^l^K-^l^ ; 

this is to be considered to be stationary in time so that it also represents the ide- 
alized time average. These bounds and the idealized energy spectral function are 
illustrated in Figure lTT] The first constraint implied by (13. 61 ) (13. 71 ) and (13.81 ) is that a 
spectral regime with parameter e is incompatible with the situation in which Ek{k) 
lies entirely above the permitted set S:={E < 47iRl} n{E < 4kRI{T)/vk^T}. 




ln(E(K)) 



Figure 3.1. The accessible set S and a spectral function Ek{k). 



Proposition 3.3. In order that the graph ofEK{K) intersect the set S, the parame- 
ters must satisfy the relation 

5/3 

(3.10) 



vV6co£2/3 < 4;r r^") r\I\t). 



The proof is elementary. This gives an upper bound on the parameter e, in 
fact on the quantity CoV^/^e^/^, in terms of quantities that are determined by the 
initial data and the inhomogeneous forces. In the setting of statistically stationary 
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solutions, R\{T)/T < R2, where R2 and R\ are constant. In order that a spectral 
regime is achieved, the relation 

must hold. This constrains the values of the parameter e for any solution regime 
that exhibits spectral behavior. 

We now take up the question of the endpoints of the inertial range [k"i , K2], as- 
suming a given value of £. We will produce an interval [Tci ,7C2] such that upper and 
lower limits of the inertial range, respectively Ki and K2 must necessarily satisfy 
^1 < K'l < K'2 < ^2- First of all, the function Ek{k) will violate the estimate ( 13.61 ) 
(if the force is not present) or (13.71 ) (when there is a force) unless k>Ki, where 

(3.11) Cae'^I^K^^'^ =AnR\, 

which gives a bound from below for the lower endpoint of the inertial range. 

Proposition 3.4, An absolute lower bound for the inertial range is given by 



(3.12) Tci- "^0 



(47r/?2)3/5 • 

It is an amusing exercise to check that the RHS has the appropriate units of 
L^', for which we note that the units of e are l} /T^. In the case of a nonzero 
forcing, R\ (t) may be increasing, in which case Tci (f ) would decrease. In the case 
of bounded forces, Ri{t) may increase linearly in t, implying that Ki{t) r-^ t^^'^. 
For a statistically stationary solution as described above, R and Ri, and therefore 
Tci are constant. 

The upper bound for the inertial range comes from comparing time averages of 
Ek with the upper bound (13.81 ). Indeed, 



(3.13) i rE,{K)dt=Coe'^'K-'/' < ^^^ . 

T Jo vTk"^ 

Proposition 3.5. The inequality (13.131 ) holds only over an interval of K bounded 
above by 

{Anf 1 R\{T) 



(3.14) K2 



(C()V)3£2 73 



It is again amusing to check that the RHS has units of L^ , noting that v has 
units of iJ- jT . When there is a force present, the constant R^(J^ may grow in 
r. When / G L~(M+;//"^ PlL^) the constant RiiJ) grows at most linearly in 
r. When considering the case of a bounded and statistically stationary forcing 
term, for example, the ratio R\{T)/T is expected to have a limit as T grows large, 

\miT^+ooR\{T) /T = R2, which gives rise to a fixed upper bound for 7c2. Indeed in 
any case in which the constant /?2 := hmsup7-^^^(/?2(r)/r) is finite, this argu- 
ment gives an upper bound for K2. 
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However with no force present, or with a force which decays in time, then 
R\{T) will be bounded, or may grow sublinearly, which results in the bound for 
7C2 = ^2(7") which is decreasing in time. Supposing that at some time Tq we have 
that for T >Tq then lC2(r) < Tci, implying that the interval consisting of the inertial 
range is necessarily empty. The explicit bound for Tq in the case of no force present 
is as follows. 

Proposition 3.6. Suppose that the force f = 0, so that Ri and R2 are constant in 
time. Then K2{T) < Kifor all T > Tq, where 

The RHS has units of time. If there is a nonzero force present, then R\ = Ri{T) 
and R2 = R2{T), so that the expressions (13.121 1 (13.141 ) for "^i = Ki{T) and K2 = 
K2{T) depend on time. It nonetheless could happen that 

(3.16) Iimsup7c2(r) < liminf7ci(r) , 

then again there is a maximum time Tq for the existence of spectral behavior of 
solutions. 

The above three estimates give lower and upper bounds on the inertial range, 
and an upper bound of the time of validity of a spectral description of a solution 
to (12. lb . if indeed it behaved exactly like the Kolmogorov power spectrum profile 
over its inertial range. 

As discussed in the introduction, when additional physical assumptions are 
made as to the behavior of a solution of (12. lb . then further information is avail- 
able about the energy transfer rate e. The classical hypotheses are stated in |[T4l 
among other places, describing the character of flows in a regime of fully devel- 
oped turbulence. Specifically, they are that (1) the flow is in a (statistically) steady 
state of energy transfer from the inertial range \k\ < Ky to the dissipative range 
1^1 > Kv', (2) the support of the spectrum E{K,t) Ues essentially in inertial range; 
and (3) a certain scale invariant form is assumed for the transport of energy T{k) 
at wavenumber scale K which assumes a form of homogeneity of the flow. Under 
hypotheses (1)(2) and (3), one concludes that the energy dissipation rate £1 in (11.61 ) 
is equal to the energy transfer rate e, and that the upper end of the inertial range 
K"2 = K"v = 2n{e/v^y'^. Alternatively, one can simply work under the hypothesis 
that the energy spectral function depends only upon the two quantities k and the 
energy dissipation rate £1 . Assuming that the force is stationary, and that £ = £1 , 
the conclusion of lH is that 

where (•) denotes either time or ensemble averaging, and in any case ||m(-)||l2 ^ 
7?. This is to be compared with the general estimate (13.101 ). and it behaves better 
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for small v. On the other hand, the estimates in the present paper hold for all 
weak solutions of (12.11 ). and they give an upper bound on the energy transfer rate 
£, essentially independently of the Obukhov hypotheses. Furthermore, the upper 
and lower bounds on the inertial range are conclusions of the analysis rather than 
assumptions of the theory. In particular this work gives a lower bound on the 
inertial range, which is the first such known, either with mathematically rigorous 
arguments or under physical hypotheses, at least to the authors. 

3.4 Limits on spectral behavior 

The endpoints K^ , 7C2 of the interval which bounds the inertial range, and the 
temporal upper bound Tq are given in terms of the idealized Kolmogorov spec- 
tral function Ek{k), rather than one given by an actual solution of the Navier - 
Stokes equations. In order to have a relevance to actual solutions, one must quan- 
tify the meaning of spectral behavior of a solution. This can have a number of 
interpretations, several of which we have mentioned in Section 13.11 It could be 
that we define an individual solution to have spectral behavior if its energy spectral 
function E{kj)i& sufficiently close to the idealized Kolmogorov spectral function 
Ek{k), uniformly over a time period [0,T]. Since the solution u{-,t) G L^ and is 
indeed in H^ for almost all times t, while Ek is neither (i.e. neither Ek{k:) nor 
k^Ek{k:) G L} (MJ)), this already implies that the inertial range must be finite. 

Definition 3.7. A solution {u,p) to (12.11 ) is said to have the spectral behavior of 
Ek{k), uniformly over the range [k"i , K2] and for the time interval [0, T] if its energy 
spectral function E{K,t) satisfies 

(3.17) sup _{l + K^/^)\E{K,t)-EK{K)\<Ci«Coe^/\ 
Ke[K[,K2],te[0j] 

An alternate version of this specification would be to replace the criterion (13.171 ) 
with a weaker one, for instance asking that a Sobolev space norm be controlled, 
which for D = M^ could be the statement that 

(3.18) sup / {l + K^'^)\E{K,t)-EKiK)\dK<Ci«CQe'^^\ 
re[Oj]"'K'i 

Or else one could specify a criterion which respected the metric of a Besov space. 

For example, one could use the Fourier decomposition Aj = {k : l^^^l"^ < \k\ < 

2^+V2|^ and ask that over a time period [0,T] a solution satisfy 

(3.19) 

/ |M(A:,OP^/t-Coe^/^3sinh((ln2)/3)2-2i/3 <ca'^il^ ^Cae^'^l-^^l^ 
Jaj 

for all 7i < 7 < J2, where j'l , 72 are such that 71 < log2 K"i , and logj K2 < ji- In any 
of these cases. Theorems 13. 1 1 and 13.21 imply bounds on the inertial range given by 
the interval [^i,7c2] . 
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Theorem 3.8. Suppose that an individual solution {u{x,t),p{x,t)) is such that 
uo{x) G As, 06^(0), where R and Ri satisfy (12.91 ) and if a force is present, it satis- 
fies (I2.12I ). Ifu{x,t) exhibits the spectral behavior ofEK uniformly over the range 

[ki,K2] X [0,T], then 

(3.20) "JCl < JCi , K2<K2, 

and iff = then 

(3.21) T<To, 

with possibly different constants Co in (I3.12I ) (I3.14I ) (I3.15I ). 

Thus the spectral behavior of solutions whose initial data uq{x) lie in one of the 
sets Ar^ nB/;(0) is limited by the bounds given in (13.201) . The proof will show that 
the same constraints hold for solutions which exhibit spectral behavior over a given 
nonzero proportion of the measure of the time interval [0,T]. 

However it could be argued that the behavior of an individual solution is less 
important, and that spectral behavior is a property of a statistical ensemble of so- 
lutions. Members of this ensemble should have their spectral behavior considered 
in terms of the ensemble average, rather than individually as above. Theorems 13.11 
and 13.21 are relevant to this situation as well. Suppose there were a probability 
measure P defined on a statistical ensemble Q. C L^{D) n {m : V • m = 0} which 
is invariant under the solution map of the Navier - Stokes equations, however this 
has been chosen to be defined, with force / (also possibly stationary, taken from a 
family of realizations which have their own statistics). Using the standard notation, 
define the ensemble average of a functional F{u) defined and P-measurable on D. 
by {F{u)). Without loss of generality we can take P to be ergodic with respect to 
the Navier - Stokes solution map. 

The ergodicity of the invariant measure P tells us two things. The first is that 
space averages are a.e. time averages, so that 

(3.22) (EiK,-) =lim- EiK,t)dt <—^lim^ 

for P-a.e. initial data uq. The second thing is that whenever R,Ri satisfy (12.121) 
then P(As| nB«(0)) is either zero or one, as Ar^ nBs(O) is an invariant set. 

Definition 3.9. A statistical ensemble (Q, P) is said to exhibit the spectral behav- 
ior of Ek{k) on average over the range [k"i , K2] when the ensemble average of its 
energy spectral function, 

(3.23) {E{K,t)):= [ {\u{k,t)\^)dS{k) 

J\k\ = K 

satisfies the estimate 

(3.24) sup _{l + K^/^)\{E{K,t))-EK{K)\<Ci«Coe^/^ 
Ke[Ki,K2],te[0,T] 

over the range [kti , K2]. 
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Let us suppose that the force / satisfies 

(3.25) |/>,0|<v|/:|/?i and {j\f{k,s)\^dsY'^<F^{t), 

as in (12.121 ). and we are to examine the spectral behavior of the statistical ensemble 
of solutions {m(-)}. 

Theorem 3.10. Suppose that the ensemble {D., P) has the spectral behavior of 
Ek{-) over the range [k"i, K2]. Then either 

(3.26) P(A/,,nB«(0))=0 
for all R,Ri, or else 

(3.27) "JCi < K-i , K2<K2 , 

with a possibly different constant Cq in ( |3.12| )( [3TT4l )( |3.15l ). 

Proof of Theorems \3. 8\ and \3. 1 0\ We will give the argument in the case of Euclid- 
ian space D = M^, the torus case is similar. Start with the proof of Theorem l3.8l with 
the criterion of (13.171 ). and suppose that K\ < TCj. Using the estimate of Theorem 
Oand the form ^^ ofEK{Ki) we have 

Coe^/^i"^^^ -47r/?? < \Ek{ki)-E{ki)\ < o(l)Coe^/^ . 
Because of the identity ( 13.111) . this implies 

Coe'/\K;'/'-K-'/')<o{l)Coe'/'. 
A lower bound for the LHS is given by 



-5/3 --5/3 



£lK-'/'dK>Oc,-K,)lK-'^\ 



)-7c'/' 



Therefore 

0<(7ci-fCi)<o(l)|7c° 

that is, K"i is a bounded distance from Tci. Furthermore the defining relation (13.111 ) 
for the left endpoint "^ of the bounds on the inertial range can be rewritten 

4nRl = Coey%'^' =C,ey'K;'/'[^f , 

and since (1 — o(l)(3/5)7c/ ) < Ki/K[ < 1, then ( 13.111 ) continues to hold for Ki, 
with only a small change in the constant Cq. 

Now suppose that K2 ^ ^2- The criterion (13.171 ) implies that 

U\'^'\EK{K)-E{K,t)\dt<o{\)Coe^/\ 
T Jo 



Therefore using (13.81 ) and ( 13.91) . we have 



Coe2/3_l!EM^-l/3<o(l)Q,£2/3 



26 A. BIRYUK & W. CRAIG 

for K" G [k"i , K2\ and T <T. This applies in particular to K" = K2, therefore 



V T ^1/3 ■ 
Hence 



^2 



where we have used ( I3.14I ). Therefore 

(l-o(l))<^<l, 

and (13.141) holds for K'2 with only a change of the overall constant Cq. 

Similar considerations give the analog result to Theorem 13.81 if we accept the 
Sobolev or Besov criteria for spectral behavior. For instance, suppose it is con- 
sidered that the estimate (13.181) is the indicator of spectral behavior. If Ki < Ki 
then 

/ \coe^/^K-^^^ -47iRJ)dK < [ '\EK{K)-E{K)\dK<o{l)Co£^/^ . 

Using (13.11b . this implies that 



Jki 



,2/3 

which in turn implies (by convexity) that 



This controls Tci — K^i and also their ratio. Suppose that K2 < K2. Then the criterion 
(13.181) impUes that 

r K'/\Coe^^'K-'/'-^K-^)dK< r K'l^\EK{K)-E{K)\dK<o{\)Coe^/^ 

Jki VT J K2 

Therefore, using ( 13. 14b 

(k-2 -K2) - -K^^^K^^^ -kI^^) < 0(1) . 

Define a{K) := k:— ^7^2 k^^^, which is increasing and convex for k>K2- The 
last estimate states that a{K2) — a(7c2) < o(l), which controls the quantity K2 — "^2- 
The proof of the analogous statements of Theorem l3.10l are similar, except that 
the upper bounds on K2 are easier as the ensemble average already subsumes the 
time average due to the ergodicity hypothesis. D 
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4 Conclusions 

The global estimates given in theorems |2.3[|2.4| and |2.5| for the domain D = T^, 
and theorems 12.71 and |2^ in the case D = M?, provide control in L°° of the Fourier 
transform of weak solutions of the Navier - Stokes equations. These are in terms 
of constants R, Ri and /?2 which depend only upon the initial data and the inhomo- 
geneous forces. These results in turn give estimates of the energy spectral function, 
which show that E{K,t) is bounded from above, and its time averages are bounded 
above by ^(1/k"^). These upper bounds constrain the ability for a weak solution 
to exhibit spectral behavior in the manner of the idealized Kohnogorov spectral 
function Ek{k) = CqE^'^k^^'^ . The constraints extend to the case of a statistical 
ensemble forces and solutions, applying to the ensemble averages {E{K,t)) of the 
energy spectral function. We remark that the estimates, and the subsequent con- 
straints on spectral behavior, are valid for weak solutions of the Navier - Stokes 
equations, and our considerations are separate from the physical assumptions of 
Obukhov on flows exhibiting fully developed turbulence, or the question of possi- 
ble formation of singularities. 

It is natural to compare the above constraints with the physical quantities de- 
scribing spectral behavior and the inertial range that come from the Kolmogorov - 
Obukhov theory of turbulence. The first of these is the Kolmogorov length scale 
rjv = {v^ /e)^'^, or rather its associated wavenumber Ky = 2n/r\y. On physical 
grounds, dissipation is expected to dominate the behavior of £(«:,?) for K" > Ky. 
Comparing Ky to our upper bounds on the inertial range, we find that 

for sufficiently small £ and V. Indeed, with everything else fixed, Ky is decreasing 
in e while 7c2 is increasing, and furthermore while both Ky and K2 are increasing as 
V — > 0, however Ky '^lC2- It seems clear that K2 is an absolute, but not necessarily 
a very sharp, estimate of the upper limit of the inertial range and the start of the 
dissipative regime for solutions that is expected on physical grounds. 

As described in section 13.31 if one assumes a number of physical hypotheses, 
as Obukhov does, on the form of the energy transfer rate, from which one deduces 
that Ky gives an upper bound on the spectral regime and that £ = £1 is given by the 
energy dissipation rate, then there are better upper bounds available for £ than our 
result ( 13.101 ). This assumption however is based on physical assumptions on the 
character of solutions of the Navier - Stokes equations in a statistically stationary 
regime. 

The Taylor length scale Kx = 2n{£V /vR^Y'^ is another indicator of the lower 
limit of the dissipative regime, one which incidentally is independent of the form 
of the Kolmogorov idealized energy spectral function Ek- The quantity R^ /V is a 
bound on the energy per unit volume of the solution. We again see that K2 is an 
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overly pessimistic upper bound for K^ for small e and v, since in such a case 

vR^J - \CoV T ) s^ 



'^A=27r(::^) < ( ^^^A_^ ) - = k2 



In both of these comparisons the quantities R\{T)/T are to be replaced by Rj in 
the case of a statistically stationary ensemble of solutions. 

In a flow regime of fully developed turbulence, it is generally expected that 
Ki <Kv, an inequality which is worthwhile to discuss. Calculate 

,3x1/4 



^X _VV/£V/Vv'A ' _ VV,^, 1/4 



Ky R \VJ \ £ J R 

where {evY^'^ = Uy is the Kolmogorov velocity scale. For solutions that we con- 
sider, (12.91 ) holds, so that in particular 

VV/_..^l/4^ ^ /_..n1/4 _ R Ky 



R ^ ' - vRi ' Ri2n 



The implication is that 



R 2 




which indicates, for fixed data R,Ri, that K^ cannot be too much smaller than Ky, 
and is very possibly much larger. An inequality in the other sense does not seem 
to arise from this or similar considerations. We do find an upper bound for the 
Kolmogorov velocity scale 

"^"^'"^^ -VvKy-vy^^[cJ [[vf 

Considering the case / = 0, which is as in the original papers of Kolmogorov fS"], 
the constraint of Theorem l3.8l is that T < To, with the latter given by the expression 
in (13.151) . This is to be compared with the Kolmogorov timescale Ty = (v/e)^/^. It 
is clear, for Ri and /?2 fixed constants, that 

^ ^ £3/10^3/2 [ ^0 I 

To \{47i)(>/5r]^'rjJ 

which is of course small for small e, v. This is again as it should be, allowing large 
multiples of the eddy turnover time before one runs into the upper allowed limit 
for the persistence of spectral behavior of solutions. 

It is also natural, given the bounds Ki , K2 on the inertial range, to introduce the 
dimensionless parameter 

K2_ 1 f4nY'/'(RfRl{T)y 



^'■'^ '^--=¥,=^^^^[-^0) 
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which governs the extent of the possibihty of spectral behavior of solutions. It is 
somewhat similar to a Reynold's number; when ry < I then solutions satisfying 
( |2.9I ) (I2.12I) (respectively, ( I2.19I )) are disallowed from exhibiting spectral behavior. 
For ry > 1 an inertial range is permitted, although it is not guaranteed by the anal- 
ysis of this paper. The larger ry the larger the permitted inertial range, although 
again it is not the case that the actual interval of k over which solutions exhibit 
spectral behavior will necessarily extend through a significant proportion of the in- 
terval Ki,K2- In the situation of a statistical ensemble of forces and solutions, the 
form of ry is somewhat more compelling, 

^'•'^ '"-= ^ = iw^;^ UJ r^ ""') • 

This quantity is a stand-in for the ratio of the integral scale to the Kolmogorov 
scale, which is itself often used as an indicator of the Reynolds number of a flow. 

This paper does not address the corrections to the Kolmogorov - Obukhov the- 
ory of Navier - Stokes flows in a turbulent regime, along the lines proposed in 
Kolmogorov (1962) 1 10|. This is focused on the deviations from Gaussian nature 
of the moments of the structure function for such flows, and it has been a very 
active area of research over the past decades. We will reserve our own thoughts on 
this matter for a future publication. 
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